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Abstract 

Transport theory for the Nambu-Jona-Lasinio (NJL) model is exam¬ 
ined. The collision term is investigated in a first approach in a coupling 
strength expansion. It can be demonstrated that this leads to a form 
easily recognized from the Boltzmann equation if the quasiparticle ap¬ 
proximation is used. It is seen that enforcing a quasiparticle approx¬ 
imation suppresses the three body creation and annihilation processes 
that would otherwise be present. Including mesons and hadronization 
consistently in a 1/Ac expansion is discussed briefly. Some comments 
on the numerical simulation of the Vlasov equation and the relaxation 
time approximation to the collision term are made. 

1 Introduction 

It is a well accepted fact that the property of chiral symmetry is essential for 
a correct description of the static properties of the low lying mesons. Fnrther, 
it is known from lattice simnlations that, at a hnite valne of the temperatnre, 
this symmetry is restored. This, in tnrn, may lead to snbstantial changes 
in the mesonic spectrnm, something that we know only from model calcu¬ 
lations |l|. Now, with the increasing interest in heavy-ion collisions, and in 
the non-equilibrium descriptions thereof, this particular feature of quantum 
chromodynamics (QGD) has to date been ignored. Thus, the purpose of this 
work is to construct a transport theory that is based on a model of QGD 
that contains the essential feature of chiral symmetry breaking and restora¬ 
tion, and to investigate by simulation whether dynamical effects due to this 
symmetry breaking / restoration are evidenced that could have experimental 
consequences. 

Our starting point is the Nambu-Jona-Lasinio (NJL) model, that intrinsi¬ 
cally contains only quark degrees of freedom that interact via an effective two 


and three body interaction for three flavors. Mesonic states are constrncted 
as collective excitations in this model. No explicit glnonic degrees of freedom 
are inclnded. A transport theory for this Lagrange density, at the lowest level, 
involves only the qnark degrees of freedom. One is reqnired to go to higher 
orders - in this strong conpling theory, an expansion in the inverse nnmber of 
colors 1/Nc is appropriate - in order to inclnde the mesonic degrees of free¬ 
dom. Doing so, one hnds that it becomes possible to constrnct a consistent 
theory of qnarks and mesons in which the collision term contains inter alia the 
hadronization processes gg' — MM' in addition to the usnal elastic scattering 
qq' —*■ qq' and qq' qq' in the quark sector. In this talk, we sketch the de¬ 
velopment of the collision term hrst expanded in the interaction strength, and 
indicate how mesons are included. Enforcing a quasiparticle approximation 
is seen to suppress the three body creation and annihilation processes that 
would otherwise be present. The expansion shows how meson dynamics are 
incorporated naturally into the one-body transport equation for the quark or 
antiquark distribution function. 

Finally, we indicate briefly the current status of the numerical simulations. 

2 Including the collision term in transport theory 

2.1 Equations of motion 

Non-equilibrium phenomena are completely desccribed via the Schwinger - 
Keldysh formalism for Green functions 0. Many conventions exist in the 
literature. For our purposes, it is simplest to use the convention of Landau [§ . 
In this, the designations -|- and — are attributed to the closed time path that 
is shown in Fig.l, and the fermionic Green functions are dehned as 
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In a standard fashion, one constructs equations of motion for the matrix of 
Green functions and one moves to relative and centre of mass variables, u = 
X — y and X = [x + y)/2. A Fourier transform with respect to the relative 
coordinate, or Wigner transform, 

S(X,p) = j d\e‘<-'‘S (v + |. V - I) 


(2) 


is then performed. Of particular interest is the equation of motion for S = 
S^. Regarding this, together with the equation of motion for the adjoint 
function one arrives at the so-called transport and constraint equations 
by adding and subtracting these. One hnds 

y { 7 ", + P,[r, S<] = ^<S> - E>S< + [S^ ^>] + [S^, S>], (3) 

and 
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which are general equations that are generic for any fermionic theory [Q. 
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Fig. 1. Closed time path on which Green functions are defined. 


2.2 Methods of solution 

A complete method of solution of Eqs.(H) and (Q) follows on inserting the 
spinor decomposition 

= F + i-f^P + YVf, + (5) 

for into the above equation. Unfortunately, this is only practical for the 
collisionless system 0, and it leads to a set of 16 x 2 equations that must be 
solved simultaneously. In order to understand the collision term, it is more 
useful to insert the Ansatz 

S^{X,p) = 27ri^[5{po-Ep)J2us'ip)usip)ffip,X) 

s,s' 

+ 5{j)q + Ep)^Vs'{-p)vs{-p)ff {-p,X)], ( 6 ) 

ss' 

with / = ! — /, and evaluate the determining equation for the quark and 
antiquark distribution functions fg and fq (note that = Sss'fq)- 








2.3 The collision term 

It is useful to understand the simple collision graphs that are given in Fig. 2, 
where the interaction is displayed as extended for convenience. 



Fig. 2. Lowest order contributions to the collision term. 


Then the terms on the right hand side of Eq. (|^) that describe collisions are 


C.I. = S>(X,p)5<(X,p) - S<(X,p)5>(X,p), (7) 


and which represent loss and gain terms respectively. For the simplest form of 
the NJL model p, three possible forms for S exist, Eg-, and Tjmixed, which 
contain, for the spinor parts of the interaction vertices, unity, fys throughout, 
or *75 and unity in a mixed fashion |p. For example. 


S?(V,p) 
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Now the insertion of Ansatz (^) can be used. One hnds that 

• the combination of all three types of self-energies is required in order to 
construct the complete relativistic scattering matrix elements for quark- 
quark and quark-antiquark scattering, and 

• processes that lead to the production/destruction of quark-antiquark 
pairs are present even at this level 0, but are zero if on-shellness is 
stipulated. 

The loss term, obtained as 
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where A+ is a restriction over positive energies is then fonnd to be 
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where M{p2 —>• 13) = M°'(p2 —13) + M'^{p2 13). Thus one has the 

relativistic generalization of the nonrelativistic calculation of Kadanoff and 
Baym |0]. At this level, only elastic scattering is incorporated. 


2.4 Higher orders in 1/A"c and meson prodnction 


Although the previous calculation is informative, an expansion in the coupling 
strength is inadmissable since GA^ ~ 2. We thus briefly delineate the extension 
to higher order, leaving the details to a forthcoming publication p|. To do this, 
we note that the expansion in 1/Nc must be performed in such a way as to be 
symmetry conserving. The next order terms have been investigated in detail 
in Ref. IP . The self-energy, in the next order, contains meson exchange, leading 
to a self-energy of the form 


m 


/ - Cl)}, 


where the superscript E describes the self-consistently calculated Green func¬ 
tion at this level, see Fig. 3. 



Fig. 3. Self-consistent self energy that includes meson exchange. 


Expanding these in terms of the Hartree propagator for S', and the lowest order 
pion irreducible polarization 

(g)-^ = 1 - 2GYf{p) = 1 - 2G(n° + OT), (12) 

see Fig. 4, enables one to obtain inter alia the diagrams for S that are shown 

in Fig. 5. 
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Fig. 4. Contributions to the irreducible polarization to next order in l/Nc- 


Intuitively it is obvious that cutting rules applied to these graphs (the hrst one, 
we term “rainbow”) lead to the hadronization cross sections as were calculated 
in Ref.|]^ and which now form part of the collision integral. Thus they are 
incorporated naturally into the transport formalism and are in fact essential 
in a consistent l/N^ expansion. Details of this calculation will appear in a 
forthcoming publication |^. 



+ 



Fig. 5. Diagrams occurring in S that when cut, lead to hadronization of a 
quark and antiquark into two mesons. 

3 Numerical simulations 

Assuming a radial symmetry, the Vlasov equation for the NJL model i.iig 


has been directly integrated using a hnite difference method [111. Starting with 


a temperature of T = 240MeV, and /i = 200MeV and a factorizable prohle 
/(dp) = frir)fp{p) with fr{r) = (1 - tanh((r - ro)/6r))/2, tq = 3.1fm and 
Sr = Ifm and fp{p) = 2NcNf /— p) + 1), we hnd that the symmetry 
breaking order parameter is globally almost constant at times not exceeding 
17fm/c, setting a maximal time for the symmetry breaking in this system. 

Calculations that include the collision integral via a relaxation time ap¬ 
proximation are at present underway. 
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